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abstract 

We study how the initial inhomogeneities of the spatial curvature 
affect the onset of inflation in the closed universe. We consider a 
cosmological model which contains a radiation and a cosmological 
constant. In order to treat the inhomogeneities in the closed uni- 
verse, we improve the long wavelength approximation such that the 
non-small spatial curvature is tractable in the lowest order. Using 
the improved scheme, we show how large inhomogeneities of the 
spatial curvature prevent the occurrence of inflation. 



I. INTRODUCTION 

Inflationary universe model Q is the most favorite scenario which can explain the large scale homogeneity and 
isotropy of the Universe. One of the main questions concerning the inflation is whether the occurrence of inflation is 
generic or not. 

There are lots of investigations on this problem in the homogeneous model universes. For Friedmann-Robertson- 
Walker universe with an inflaton ^|J|] , the occurrence of inflation is generic except the closed universe case. In the view 
point of cosmological no-hair conjecture, Wald Q showed that initially expanding anisotropic Bianchi type universes 
with a positive cosmological constant except the type IX approach de Sitter universe. In theses cases, inflation is not 
completely generic because the positive spatial curvature can prevent inflation. 

It is obvious that the most appealing feature of the inflationary cosmologies is smoothing out the initial inhomo- 
geneities of the universe. Thus, we should clarify the generality of inflation in the initially inhomogeneous universe. 
The effect of inhomogeneities on the inflation is studied by numerical simulations |||| ■ I n ^ ei ■ @ j it is shown that 
the homogeneity over a few initial horizon is needed in order to inflate the universe and other works support it. 

There are two aspects of the inhomogeneities. One is the inhomogeneities of the inflaton, the other is the one in the 
metric. In general these are coupled. In this paper we study how the inhomogeneities of the spatial geometry affect 
the occurrence of inflation by the use of the long wavelength approximation, which reduces the Einstein equations to 
several sets of ordinary differential equations. 

By the use of the long wavelength approximation method we can construct approximate solutions of Einstein equa- 
tions which describe an inhomogeneous universe on scales larger than the local Hubble radius. Belinski, Khalatnikov 
and Lifshitz Q| studied the general behaviors of the universe in the neighborhood of initial singularity by this ap- 
proximation, and Tomita || focused attentions on the evolutions of inhomogeneities of the early universe in this 
framework. 

The long wavelength approximation assumes that the spatial derivative of a metric is smaller than the time deriva- 
tive. The spatial curvature which consists of the spatial derivative of the second order is neglected in the lowest order 
of this approximation and it is taken into account perturbatively by the series expansion in the higher order. This 
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approximation is called the gradient expansion method for the Einstein equations. For the higher order solutions, 
Comer, Deruelle, Langlois et al Q developed an iteration scheme of Einstein equations and Salopek, Stewart et al 
developed one of the Hamilton- Jacobi equation for general relativity. 
Recently, using the gradient expansion method, the generality of inflation is discussed by Deruelle and Goldwirth 
JTT| ], and Nambu and Taruya Il2]. As shown by the investigations of the homogeneous models [pP], the positive 
spatial curvature plays an essential role; if the positive curvature dominates over the vacuum energy, which drives 
inflation, the universe never enter the inflationary era and recollapse to the singularity. In the marginal solution 
between inflationary and recollapsing solution, the spatial curvature becomes as large as the vacuum energy. For 
desired inflation to explain the homogeneity of the universe, large vacuum energy is necessary. Then we should take 
account of non-small spatial curvature to clarify the condition for the occurrence of inflation. Readers would suspect 
that the gradient expansion method is not suitable for investigation of the generality of inflation because the spatial 
curvature is assumed to be a small quantity. 

In this paper, we improve the gradient expansion scheme so that the non-small constant spatial curvature can be 
tractable in the lowest order. By this improved gradient expansion scheme, we study how the inhomogeneities of the 
spatial curvature affect the occurrence of inflation in the closed universe. The word 'closed universe' does not mean 
the universe with spatial surface of closed topology in this paper. We investigate the time evolution of the spatial 
region with non-small positive spatial curvature. 

For the first step to investigate the effect of the inhomogeneities of the spatial curvature for the onset of inflation we 
consider the model in which the cosmological constant gives rise to inflation instead of the scalar field and a radiation 
fluid fills the universe. 

This paper is organized as follows. In Sec. 2, we explain the improved gradient expansion method in order to 
construct the inhomogeneous closed universe. And we derive the basic equations in the first order approximation. In 
Sec. 3, we specify the basic equations for our simple model and show numerical results for our aim. Sec. 4 is devoted 
to conclusion and discussions. 



II. GRADIENT EXPANSION FROM CLOSED FRIEDMANN UNIVERSE 



In this section, we review the usual gradient expansion method briefly and modify it so that non-vanishing spatial 
curvature is included in the lowest order. 



A. Brief Review 



We start from a brief review of the usual gradient expansion method. 
Take a synchronous reference frame where the line element is 

ds 2 = -dt 2 + jij(t, x k )dx i dx j . 

Throughout this paper Latin letters will denote spatial indices and Greek letters spacetime indices. The matter is 
taken to be a perfect fluid characterized by energy-momentum tensor 

Tftv = {p + p)u^u u +P9nv, 

where p, p and u M are pressure, energy density and four velocity, respectively. We assume the equation of state of the 
fluid is p/p = r — 1, where T is a constant. The Einstein equations with a cosmological constant A are 

2K + KTKL-4A 

KP ~ 2(2 - 3r - 2Tu,u') ' [Z - Z) 

nTpui = -—L=(Ki tj - K i), (2.3) 
2V1 + W 

where 3 i?^ is the Ricci tensor associated with 7^, Kij = 7^ is the extrinsic curvature, 7 = det7y, a dot denotes the 
derivative with respect to t, a semicolon denotes the covariant derivative with respect to 7^ and k = %ixG. 
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We consider the situation in which all spatial gradients are much smaller than time derivatives and assume that 
the spatial metric is expanded as 



(o) (i) 
7y (t, x)=l ij (t, X)+ 7 ij (i, x) + 



lij (t,x) 

p 



(2.4) 



(0) 



Here the lowest term 7 y in Eq.(2.4) is assumed in the quasi isotropic form 



7 ij (t, x) = a 2 (t)h i: j{x), 



(2.5) 



where a(t) is the scale factor in the lowest order which depends on only t and hij{x) is an arbitrary function of only 

w . 

spatial variables. The next-order term 7 ij is constructed by the symmetric tensors of the second order in the spatial 
gradients of h{j in the form 



7ij (t,x) = a 2 {t) [/i(t)i2y(fc) + g x {t)R{h)h 



ij J ■ 



(2.6) 



where Rij(h) is the Ricci tensor associated with hki(x), R(h) = Rij(h), and fi(t) and gi(t) are functions of time. 

( P ) 

The subsequent general terms 7 ij in Eq. fl2.4| ) are spatial tensors which are of order 2p in the spatial gradients of 
hijix) with time dependent coefficients. Since the spatial inhomogeneities of the spacetime is generated by hij(x), 
then we call it 'seed metric'. 

Substituting Eq.([2.5|) into Eqs.(2.1)-(2.3) and neglecting the spatial derivatives we get the lowest (zeroth) order 
solution 



a(t) = 



sinh 



3r 



-H Q t 



(o) 
P : 



3H§ 



Ksmh 2 (^-H t) 



(0) n 
u i= 0, 



(2.7) 
(2.8) 
(2.9) 



where Hq = A/3. The evolution of scale factor at zeroth orde r is th e same as the one of flat Friedmann-Robertson- 
Walker universe. So, we call the solution described by Eqs.(2.7)-( |2.9| ) "locally flat Friedmann solution". 

(o) (l) 

At next (first) order, substituting the spatial metric 7^ =1 ij + 7 y given by Eqs. fl2.5p and ( |2.6| ) into Eq.(pj|) 
and taking the spatial derivative terms of the seed metric, we get ordinary differential equations for the functions of 
time fi(t) and <7i(t)(see ref. fl). S imilarly, we can calculate the higher orde r so lutions iteratively. The locally flat 
Friedmann universe, Eqs.(Ej)-d2j), is used as a starting point to solve Eq. ( |2 . l| j in the expansion scheme, thus we 
call this scheme the gradient expansion from the flat Friedmann universe (GEFF). In GEFF, the spatial curvature is 
treated perturbatively as a small quantity in the series expansion. 



B. Improved method 

In this subsection we improve GEFF such that we can treat non-vanishing three curvature which is discarded at 
lowest order. We show that when the seed metric has maximally symmetry, the closed or open Friedmann solution is 
derived by summing up all terms in GEFF. After that we extend the method such that the inhomogeneities around 
the closed or open Friedmann universe are tractable. 

First, we consider a seed metric in the maxmally symmetric form with non-vanishing curvature: 

hij(x) = hfj(x), 
where the Ricci curvature associated with the metric is 

i? 4J (/ l ± )=±2/ l ±. 

Since all spatial tensors of rank two derived from hfj are proportional to hfj itself, then time dependent part in 



Eq.(2.4) is factorized in the form 
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Hj(t,x) = a 2 {t) [h% + ftWRijih*) + ff i(i)J?(/i ± )/i^ + f 2 (t)R(h ± )R lj (h 



<(t)hfAx) 



where 



i(t) = a 2 (t) 1 ± {2/i(t) + 6 5l (i)} + {12/ 2 (i) + •••} + 



Substituting Eq.(2.10) into Eqs.(2.1)-(|2.2|) we obtain the equation which a± must satisfy: 



d± 2 — 3r / a± 
2~ 



si 1 T 



2 - 3r r, 

n = "A, 

2<4 2 ' 



1 



1 



±— = -(k P + A) 



(2.10) 
(2.11) 

(2.12) 
(2.13) 



These are nothing but the equations for the closed or open Friedmann universe. Of course, the universe described 
by the solution satisfying these equations is far from the flat Friedmann universe but no more than the homogeneous 
universe. 

Next, we consider the inhomogeneities of the spatial curvature. Here we think a simple case where the inhomo- 
geneities are expressed by one spatial function. We take the seed metric to be conformal to the metric of constant 
curvature, 



h ij (x)=fi 2 (x)hf j (x), 

where Q(x) is an arbitrary spatial function. 

In this case, the Ricci and scalar curvature are given by 



(2.14) 



Riiihuix)) =±2hfjx)- 



ViVj-n v t nv j n ViV l n 



- .jv o 

1 , - . V, 



Q 2 



o 



h 



(2.15) 
(2.16) 



where V$ is the covariant derivative with respect to hf 1 -. The spatial metric, therefore, in the expansion form Eq.( [2.4| ) 
becomes 

7 tf (t, x) = a 2 (t) [hij + h(t)R t] (h) + 9l {t)R{h)h t] + f 2 {t)R(h)R ij (h) + • • •] 



a 2 (t) 



{n 2 ± (2/i(t) + 6 5 i(t)) + rr 2 (i2/ 2 (t) 



■}h 



-{(h(t) + 40x(t)) + r»- 2 (±i4/ 2 (t) + ■■•) + •■■} 

+{2. 9l (t) + fi- 2 (±4/ 2 (t) + ...) + ■■■} 



fc « 



-{/!(t) + fi- 2 (±6/ 2 (0 + •••) + •••} 



+{2/ 1 (t) + fi- 2 (±12/ 2 (t) + ...) + ...} 



fi2 



+ {fr 2 (4/ 2W + ...) + ---}pp) + 



(2.17) 



Here we consider the situation that the spatial derivative of fl is small. From Eq.(2.16), it means that the derivative 
of spatial curvature is smaller than the value of curvature itself. 
Introducing 



o|(t, x) = a 2 (t) [n 2 (x) ± {2/i(t) + 6 Sl (t)} + fi- 2 (x){12/ 2 (t) + •■•} + ■ 



(2.18) 



we rewrite Eq.(2.17) as 
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Jij(t,x) = a±(t,tt) 



h tj + E F m (*. m^^tj + E F w (*. ^)(v (4) ^) 



E^)M)( v(3f)() ) 



(2.19) 



Here the notation (V^ 2p ^fi)^ denotes symbolically symmetric spatial tensors which contain 2p spatial gradients of 
Sl(t, x), where the suffix A distinguishes the tensor belonging to the same class, and F^Jt, 51) is a function of t and 
Q(x) which is determined by the Einstein equations later. This is nothing but a re-summation of the usual gradient 
expansion. We rest rict o urselves to con sider the situation which the series expansion Eq.(2.19) converges. 

Sub stituting Eq.( 2.19 ) into Eqs.(2.1) and (2.2) and neglecting the spatial gradients of fl we obtain Eqs.( 2.12 ) and 
(2.13) again for a± (t, x) for the zeroth order solution 

%(t,x) = aiit^x^ix). (2.20) 

In contrast to Eq.(2.11), the behavior of a±(t, x), defined by Eq.(2.18), depends on Q(x). The solution of Eqs.( [2.12] ) 
and ( p,13| ) is the same form of closed (open) Friedmann universe. The new expansion scheme Eq.(2.19) is based on the 
local closed (open) Friedmann universe, so we call it "the gradient expansion from a locally closed (open) Friedmann 
universe" (GECF (GEOF)). 

Next, we consider the parts of the second order in spatial gradients of il(x) in Eq.(2.1£). The tensors (V*- 2 ^)^ are 
four types and 7 y(i, x) is in the form 



lij{t,x) = a 2 ±(t,x) 



1 



G(t, n) 



n 



(2.21) 



where ViVjfl/D, and V jf2V jfl/fl 2 are defined by 



V. 



~~ -ht:- 



n 2 



v 2 



n 2 



From Eq^^) we obtain 



(2 - 3r) 2 (a' ± /a± - a ± a' ± /a 2 ± ) Q — F 



2r [3a±/a± - A] 



n 



(2.22) 



where the prime denotes the differentiation with respect to O. 
(o) , (i) 



Substituting 7^ = 7^+ 7 * and Ui given by Eqs.( p.2C| ), fl2.2l| ) and ( p.22j ) into Eq.( |2.lD , and comparing the 
coefficients of the second order derivative terms of fl, we obtain 



a± - (2- 
F + 3T— F = - — 
a± 

a ± ■- 2 
F + 3—F = — 

a± a z 



3T) 



±F-F + 2 



a'±n 



a± 



± 



a'±Q 



2(±F-F) 



G + 3T -^G = - 



_ a± 

(4-3r)(2-3r) 



G 



a± 



a± 



2ra± (3a± — Aa± 

(2 - 3r) 



F-2 [=2.- 



a±a ± 



a± 



±G + F - F'il ± 2.F 2 



F 





2 












a" 
a ± 












a± 


4A 



(2-3r) 



Ta± (3a± — Aa±) 
2 



F-2 



f2 



« 2 



±2G + 2F - 2F'fi - (±2 + -)F 2 - 2 



a' ± n 

a± 



a'lft 2 
a± 



(2.23) 
(2.24) 

(2.25) 



(2.26) 
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where we put G = G + F 2 and G = G — \F 2 for simplicity of the equation form. To solve Eg. ( [2.23 )- Eq.(2.26), we 
need to know the time evolution of F'. Differentiating Eq.(2.24) we find that F' must satisfy the equation 



F 



«± 



a 



± 
a ± 



V. 



a ± 



-4-^(±F- 



3F 



d±a± 



F) + 2(±F' - F') 



(2.27) 



III. INFLATION IN THE INHOMOGENEOUS UNIVERSE WITH A POSITIVE CURVATURE 

In this section, we consider how the inhomogeneities of spatial curvature affect on the occurrence of inflation. For 
this purpose we think a simple model, in which the inflation is driven by a cosmological constant A. We construct 
the inhomogeneous universe by the gradient expansion from the closed Friedmann universe which is introduced in the 
previous section. The matter is taken to be a ra diatio n flui d with T = 4/3 for simplicity. 

We get the lowest order solution of Eqs.( [2.I2"| ) and ( |2.I3| ) for the system in the form 



a\(t,Q(x)) = 



(o) 

Ka5_ 



2H 2 



[2Hfa 2 smh(2-ff *) + 1 - cosh(2# i)] , 



(3.1) 
(3.2) 



Here we see the behavior of the locally closed Friedmann solution Eq.(3.f ). The solution has an early phase during 
which the local scale factor behaves as a + oc t 1 / 2 . After the period the behavior of local scale factor depends on the 
local value of Q(x). The local scale factor a+(f, x) with f2(x) > fl cr approaches de Sitter solution and a+(t,x) with 
f2(x) < fl cr recollapses to the singularity, where fi cr = y / l/(2H 2 ). The local scale factor a+(t,x) with the critical 
value of Q — fl cr approaches a static solution 



3 (o)_ A 
2A' 



(3.3) 



In the early stage t « Hq 1 , the local scale factor is proportional to ^(x)^/ 2 then the spatial curvature 3 R is 
proportional to l/(£l 2 t). The fate of local scale factor is determined by the value of 3 R(t, x)t in the lowest order. The 
condition for the occurrence of inflation is 



I I 

1^ < fK 



(3.4) 



In other words, it is 



Rinit < Rcrif) — „ 9 rT j 



(3.5) 



where 3 Rinit is the spatial curvature in the early phase t « Hq 1 . 

Next, we consider the first order corrections which are proportional to the second order of spatial gradients of Q(x) 
in the form 



7tj(t,x) = a 2 + (t,Q{x)) 



-h+ + P(t,n) 



n 



(3.6) 



From Eqs.( [2.23D -( [2.27D , the equations for F,F,G and G reduce to 



G 



F + 4 



F + 3 



a 
a 
it 
a 
a 
a 
a 

"+ 



F = 



2 

7X 



F — F 



2fl 2 sinh(2H t) 



a + i ( — ) a = 4- 



F 



2 



G + F - F'ft - 
fl 2 smh(2H t) 



fl 2 smh(2H t) 



3fl 2 sinh(2F t) 



a + >! ( — ) a = 4- 



+ 



2G 
1 



2F - 2F'0 - 



3F 2 ft 2 sinh(2F t) 



1 - 



3fl 2 smh(2H t) 



2H$fl A a e + 



2fl 2 (cosh{2H t) - 1) - H a%F 



For F' we get 



-F = 



417 sinh(2F t) 



1 - 



2fl 2 smh(2H Q t) 



3F—^-(cosh(2H Q t) - 1). 



(3.7) 
(3.8) 
(3.9) 

(3.10) 
(3.11) 



We solve Eqs.(3.7)-(3.11) by numerical integration. These coupled ordinary differential equations have growing 
and decaying solutions. In the course of time, since the growing mode dominates the solution we concentrate on the 
growing solution. Then we choose initial conditions for F, F, F', G, G in the form 



F = F' = G = G = 0. 



(3.12) 



The growing mode is generated by the source terms in the right hand side of Eqs.( |3.7[ )-(3.11). For the regularity, 
asymptotic behaviors of these variables should be 



F 



1 



n 2 H ' 



t, F 



30 2 F ( 



■t, F' 



t, G 



1 



2n 2 H, 



■t, G 



r sn 2 H 



t. 



(3.13) 



Time evolution of F and G is shown in Figs.|| and ||. In the case of ft > 17 cr , F and G grow for a while and 
become constant as the universe expands exponentially. And in the case of fl < £l cr , they grow and diverge and the 
approximation breaks down. 

We demonstrate a time evolution of the spatial scalar curvature given by Eq.(B3) (see Fig.||). When we fix initial 
time, there are infinite number of combinations of parameters (fi, V;V z f2/r2, VjfiVjfi/fi 2 ) which give the same value 
of 3 R. Time evolution of 3 R depends on the parameters. Fixing f2, we vary V^Viil/il 2 and V;V ; r2/f2 keeping 

(o) 

3 R constant. When fl is fixed, since the evolution of 3 R depends only on fl, the time evolution of 3 R depends on 
ViflVifl/fl 2 . We see from Fig.|| that when ViflVifl/fl 2 is large 3 R grows rapidly. 

By the use of the trace part of the corrections, we define an effective local scale factor which includes up to the first 
order by 



a e 0, x) = 



dct 



;<>).. to. 

7 ij+ 7 ij 



1/6 



o+(t,n) 



^F(t,fl) 



WiW l fl 



(3.14) 



The effective local scale factor a e g(t, x) describes how the universe expands at the point x. 

We can follow the evolution of the effective local scale factor a e ff(t, x) by the evolution of F(t, fl) and G(t, fl). As 
is already seen, at lowest order the behavior of local scale factor a+(t, fl) is determined by the value of fl(x). On the 
other hand, the spatial derivatives of fl affect the time evolution of a e g(t, x) at first order. It means that the effective 
local scale factor evolves under the influence of the local value of fl and the value of it in the neighborhood. If a point 
in the three dimensional parameter space (fl, V;V fi/fi, W iflW ifl / fl 2 ) is specified we know the time evolution of the 
effective local scale factor a e g(t,x). 

The effective local expansion rate and acceleration rate are given by 



a cff 
a eff 
' d eff 
a cff 



a+ 
a+ 

a+ 
a+ 



F(t,fl) 



V,\7 l fl 



V. 



G(t,fl) 



\7,fl\7 l fl 



fl 2 



F(t,fl) 



2^±F(t,fl) 
a + 



WiW l fl 



G(t,fl) 



2^±d(t,fl) 
a + 



ViflV l fl 

n 2 ' 
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When the correction terms which contain the second order derivatives grow to be as large as the lowest order term, 
the expansion scheme breaks down. Since the first order terms in Eq.( [j.lE ) contain all of the second order spatial 

(o) 

derivatives then when (d+/a+) 2 > 3 R, GECF reduces to GEFF. We, thus, assume the approximation is valid while 



det(7^+ itj) 



det ifi 



det Jfj 
3 R 



c{(a + /a + )2,3$ } 



< 0.5, 



< 0.5. 



(3.15) 
(3.16) 



We divide the three dimensional parameter space (SI, V; V O/O, V;OV;0/0 2 ) into two regions: inflationary region 
and recollapsing region (see Figs.|| - ||). The effective local scale factor with the parameters in the inflationary region 
enters the accelerating expansion phase: a> e ff/ a e ff > and a> e ff/a e ff > 0. On the other hand, the effective local 
scale factor with the parameters in the recollapsing region enters recollapsing phase: CLpff/a P ff < and a p ff/a n ff < 0. 



Between these two regions, there emerges a fuzzy region where error becomes large and Eq.(3.15) or Eq.(3.16) does 
not hold before the effective local scale factor enters the inflationary or recollapsing phase. 



IV. CONCLUSION AND DISCUSSIONS 



We studied how the inhomogeneities of spatial curvature affect the occurrence of inflation in the inhomogeneous 
universe which is constructed by the gradient expansion from the locally closed Friedmann universe, in which the 
inhomogeneities are described by a function Sl(x) and its spatial gradients. For the first step, we considered a simple 
model universe with a cosmological constant and radiation fluid. 

In this paper, we expanded the spatial metric up to the second order spatial derivatives VjVjO/O and VjOVjO/O 2 . 
We solved the ordinary differential equations for the coefficient functions of these derivative terms numerically and 
evolved the inhomogeneous spatial metric |f|] . We assumed that the inhomogeneities were generated by one function 
Sl(x) but the universe needed not to have any symmetries. 

How the inhomogeneities described by V/V'O/O and V;OV'0/0 2 affect the onset of inflation is shown in Figs.|| 
- H The term V/OV'Sl/Sl 2 tends to prevent inflation and positive (negative) ViV l Sl/Sl helps the universe to inflate 
(collapse). This is related to the fact that the spatial curvature grows rapidly when V;OV z O/0 2 is large. 

From Figs.f| - 0, we get the criterion for the occurrence of inflation in the vicinity of Sl/Sl cr — 1 = V;V'£!/0 = 
V 'iSlW l Sl/Sl 2 = in a simple time independent form: 



1 

O 2 



sT~ 



V/OV'O 

o 2 



< 



n 2 ' 



(4.1) 



where a ~ 5/4 and /3 ~ 6. If the condition Eq.( [4.l| ) hol ds, t he effective local scale factor enters the inflationary phase. 
By the use of Eq.(B3) in the limit t — ► 0, the criterion (4.1) is translated into 



Ri 



( 3 R 



init ) -i 



3 R, 



3 i? 2 



(4.2) 



where a ~ 9/ 5 and (3 ~ 4 and 3 R cr (t) is given in Eq.( |3.5| ). Equations ( |4.l| ) and (4.2) correspond to Eqs.( |3.4| ) and ( ^.5] ) 



Equation(4.2) means that the occurrence of inflation at a spatial point is determined by 3 Ri n it and its gradients. The 



spatial curvature at an early time t should be less than the critical curvature value 3 R cr (t) over a region which has 
several times size of the local curvature radius ( 3 i?imt(i)) for the onset of inflation in the future, where 3 R cr (t) is 
the marginal value for inflation of the closed Friedmann universe. 

If the homogeneous universe enters the inflationary phase, it never recollapse. However, for the inhomogeneous 
universe even if the local scale factor enters the inflationary phase, it would be possible to recollapse. It may be 
also possible that the recollapsing local scale factor re-expands again in the inhomogeneous universe. In order to 
observe the final fate of the inhomogeneous universe, we should sum up all order terms of the series expansion. It is 
actually very difficult. In this paper we investigated the be havio r of s olutio ns up to the first order of the expansion 
while the criteria of the validity of the approximation Eqs.( 3~15| ) and ( 3.16 ) hold. For the purpose to know how the 
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inhomogeneities affect the onset of inflation, it is worth clarifying the behavior of the approximate solutions of the 
order. 

In the more realistic model of the inflationary universe, a certain field variable, e.g. a scalar field, drives inflation. 
In this case, not only the spatial curvature but also the field has inhomogeneities. We will discuss this case elsewhere. 
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APPENDIX A: RULE OF THE ORDER OF THE SPATIAL TENSOR 

When we calculate the first order approximation, we determine the order of spatial tensor which contains the sixth 

order terms in the spatial gradient at least because they appear in the spatial curvature 3 R/ ■ Here we explain the 
rule we adopted when we calculate the first order approximation. 

For scalar: (V (2) 0) and (V (4) 0) are the scalar of the second order terms in the spatial gradient and the fourth 
order terms, respectively. (V^O) consists of two terms, 

V;V ( V;OV'0 



o ' o 2 ' 

(V^O) is classified into the following three types, 

v,v'(v (2) o), v*(v (2) o)v'o, (v (2) o)(v (2) o). 

For spatial tensor: We adopt the tensors (V^fi)ij consist of the following terms, 

ViVjVjV'fi 

o ' 

ViVjfiVjV'ft VjfiVjVjVto VjViVj-nv'n V;V.,-(V;Ov'o) 

O 2 ' O 2 ' O 2 ' Q2 ' 

ViQVjfiVjV'n ViVj-nViftv'f} V;Ov.,-(ViOv'o) 

v i ov j ov ( ov'o 

O 3 ' 

because the trace of these terms belong to the above scalar category. 

«. 

The sixth order terms in the spatial gradient appear in the spatial curvature 3 i?/ are transformed into the following 
terms by use of the commutation relation, 

V;VmV»Vj-OV'V m O V;V t V J OV i V m V m O 

o 2 o 2 

V»Vj(V;VmV m OV i O) V t V 3 (V;V'(V m OV m O)) 

o 2 o 2 

We treat them as higher order terms. They don't affect the calculation in the first order approximation because their 
trace contain the spatial derivative of (V^O). 

APPENDIX B: EXPRESSION FOR THE SPATIAL CURVATURE IN THE FIRST ORDER 

APPROXIMATION 

In the gradient expansion from the locally closed or open Friedmann universe, we expand the spatial metric up to 
the second order of the spatial gradients in the form 
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_ (0) ±. (1) ± 



where 



(0), 



(i) 



oi(t,n(x))^, 
i 



Ifj = c4(t,fi(a:)) 



-hf i + F(t,ny 



n 



From the inverse of 7^ , we get 



It = 



4 



1 „V/V'fi,,- 



3^ 



+-(r; + 2^ 2 )^^^ + (G- J F 2 ) 



The spatial Ricci and scalar curvatures are 



1 (1) 



(1) 



3 (1) / N 



4a' ± 2^, 
3 a± 3 



a± 



a'±n 

a± 



a', SI 



+ F'fl-F + ^—F~F 
a± 



«2 



+ 



a'±Sl 
a± 

a± 



2F 



+ 2 



n 

a'±n 

a± 



+ 2F'n -2F- 



a'±fi - 5F 2 
— — F H 

a± 2 



V,OV,-0 



r> 2 



4 

2 

3 



±F-F + 2 



a± 



5? 



( 2 £ 


4\1 




4l 



ft 2 



si 



a'±Q 
a± 



±2F-2F 







+ 



±2G + 2F - 2F'Q - (±2 + -)F 2 



a'±n 

a± 



+ 



a± 



ft 2 



J 



4< 2 



±F-F + 2- 



+2 



±G + F-F , (l±2F 2 ~F 2 + n 2 2^-^ 

a± a± 



ft 2 



where a4 = da±/dCl. 
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Figure Captions 

Fig§ 

A typical example of the time evolution of F(t). 
Fig.|. 

A typical example of the time evolution of G(t). 
Figf 

Evolution of the spatial curvature which has the same initial value but different 
Figs j- 1 

An inflationary region and a recollapsing region are shown in the three dimensional parameter space 
(f2, V; V f2/0, V;£lV'f2/f2 2 ). A meshed region is a region where the approximation breaks down. 
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FIG. 2. 
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FIG. 3. 
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